Keeping a Quantum Bit Alive by Optimized 7r-Pulse Sequences 
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A general strategy to maintain the coherence of a quantum bit is proposed. The analytical result 
is derived rigorously including all memory and back-action effects. It is based on an optimized 
7r-pulse sequence for dynamic decoupling extending the Carr-Purcell-Meiboom-Gill (CPMG) cycle. 
The optimized sequence is very efficient, in particular for strong couplings to the environment. 
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Quantum information processing is a very promising 
and very challenging concept [l| . The basic feature which 
makes quantum information conceptually more powerful 
than classical information is the quantum mechanical su- 
perposition principle. It allows for the parallel processing 
of many classical registers - the so-called quantum par- 
allelism. The single quantum bit (qubit) is a two-level 
system which we may identify with a S — 1/2 system 
with states J, and |. Henceforth, we will use this spin 
language to characterize the qubit. In order for this idea 
to work the qubit has to maintain its quantum state not 
only with respect to the state f or J, but also with respect 
to its relative phase. Unavoidable couplings between the 
qubit and the environment spoil the quantum state: the 
qubit loses its coherence. This decoherence is one of the 
most serious obstacles on the way towards applications. 
Hence finding strategies to suppress decoherence is a cru- 
cial field of research [l| . 

Dynamic decoupling 0, S 0, 0| is one means to fight 
decoherence. The idea comes from spin-echo pulses in 
NMR where a large ensemble of spins is considered. 
Static, but non-uniform couplings can be compensated 
perfectly by a single 7r-pulse in the middle of the elaps- 
ing time interval. The detrimental effect of more compli- 
cated perturbations like dynamic interactions with the 
environment can be suppressed by periodic 7r-pulses or 
by periodic Carr-Purcell cycles of two 7r-pulses each Q . 

The aim of this Letter is to show that an optimized se- 
quence of 7r-pulses suppresses the decoherence even more 
efficiently than the so far known sequences of equidis- 
tant pulses 0, 0, 0, ll . The proposed scheme extends the 
known Carr-Purcell-Meiboom-Gill (CPMG) cycle [1,[3|- 
In particular for a strong coupling to the environment, 
the optimized sequence achieves a much better suppres- 
sion for the same number of pulses. So the optimized 
sequences will help to come closer to the realization of 
quantum information devices. 

We consider a fully quantum mechanical model 



(1) 



The environment is represented by a bosonic bath with 
annihilation (creation) operators b\ . The constant E 
sets the energy offset. The relevant bath properties are 
given by the spectral density d, Q 



J(w) = ^2\^\ a 8(u-Ui) 
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where we have chosen the standard ohmic bath with lin- 
early rising density in (|2b|) ; a is the dimensionless pa- 
rameter controling the coupling between qubit and bath. 
But our scheme proposed below can be applied to any 
spectral density J(u>). The high-energy cutoff ujy> is cho- 
sen as in a Debye model for phonons, but other choices 
are equally possible. Note that the correlation time tc of 
the bath is set by l/u>n- 

The model ([I]) can be easily diagonalized by the 
unitary transformation U — exp(a z K) where K — 
^ i (Ai/(2w i ))(foJ — bj) is antihermitean. The resulting 
effective Hamiltonian UK H ci = £\ u i b\b l + AE is man- 
ifestly diagonal with AE = E— J Q J(u>) /uidu). In spite of 
this simplicity ([1} suffices to study decoherence of the T^- 
type in the NMR-language which corresponds to phase 
decoherence in the XK-plane of the impurity spin. In this 
sense, ([TJ) constitutes a minimal, but fully quantum me- 
chanical, model to investigate decoherence phenomena. 
Spin flips, however, do not occur so that T\ is infinite. 
The minimal model renders the analytic examination of 
various pulse sequences possible. All thermal, quantum 
or memory effects in the bath as well as back-actions of 
the qubit on the bath are included. The pulses used in 
the following will always be considered to be ideal, i.e. 
instantaneous (cf. Refs. 0) and without any error. 

First, we look at a simple measurement assuming ini- 
tially cr 2 = 1 and the bath to be in its thermal equi- 
librium. Such a state is generated by applying a suf- 
ficiently strong magnetic field in z-dircction. Then a 
rotation about a x is applied D x (jy) — exp(—vya x /2) — 
cos(7/2) + ia x sin(7/2) which transforms the spin in z- 
direction to 



consisting of a single qubit interpreted as spin S = 1/2, 
whose operators are the Pauli matrices a x , o~ y , and a z . 



D x {^<r z D x ( 7 ) = 



cos 7<7 2 + sin 7cr y 



(3) 
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For 7 = 7r/2 a rotation by 90° is achieved; for 7 = ir the 
inversion a z — > — <r z . Measuring cr a leads to the signal 

s(t) - (t |£) 2 .(7r/2) t exp(ii?i)(Ty exp(-ii?i)£) x (7r/2)| |) 
= (T |Df(7r/2)t<(i)DS f (7r/2)| T) • (4) 

The brackets stand for the thermal expectation value 
of the bosonic bath. To obtain the line we trans- 
form by U to the effective variables and define A(t) := 
exp(iHt)Aexp(—iHt). The spin content of the resulting 
expression can be calculated using 

a c J(t)\Vl) = exp(T2Jf(t))U/T> (5a) 
af{t)\ T / |) = ±iexp( T 2K(t))\ I / T) • (5b) 

The bosonic expectation values of exponentials are com- 
puted for operators A, B linear in the bosonic vari- 
ables with the help of exp(y4) exp(S) = exp(A + 
B)exp([A,B]/2) and with the help of (exp(A)) = 
exp((>l 2 )/2). In this way, we arrive at 

s(t) = cos(2<p(t)) exp(-2x(*)) with (6a) 

**) =ifj( W )^ (6b) 

X(t) = r° JM Sm ^y 2 ^ coth(l3uj/2)du . (6c) 
Jo w 
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FIG. 1: Signal (|4j) V s. time. Solid lines for T = 0; dashed ones 
for T = O.Ic^d- Panels (a) (linear) and (b) (double logarith- 
mic) from bottom to top for a = 0.25,0.1,0.01,0.001. Panel 
(c) (double logarithmic) displays 1 — s(t) for the same values 
from top to bottom. 

Fig. Q] illustrates the effect of the coupling strength a 
and of finite temperature T = 1//3 > 0. Panels (a) and 
(b) display the usual long-time decay while (c) focusses 
on the deviation 1 — s(t) from unity. For quantum in- 
formation processes Fig. ft); shows the relevant data since 
1 — s(t) should be as low as possible. If error correction is 
to be applied thresholds between 10' " 4 and 1CT 2 

flil HU have to be met. Inspecting Fig. \Vp we conclude 
that for values of a of about 0.1 the qubit can be stored 
only for tiny fractions of the correlation time tc. But 



even if a is significantly smaller no storage is possible for 
tc, let alone for any time longer. 

Another interesting conclusion is that low values of wd 
are favorable since they set a long time scale. [23( This 
means that an elastically soft environment, for instance 
in an organic compound with low wd, is better suited 
than a hard environment with high wd- This is coun- 
terintuitive, since one might have suspected that the in- 
fluence of vibrations is lower when the spring constants 
oc are higher. The objection that the positive effect of 
a larger tc in a soft medium will be thwarted by a large 
value of the coupling a will be invalidated below. 

We pass now to a sequence of pulses where the to- 
tal time interval — > t is split into smaller intervals 
— > Sit — > S 2 t — > . . . — > 5 n t — > t with < Si < 
82 < ... < S n < 1. The S- values are taken to be fixed. 
At each instance Sit a it- pulse o~ y = —iD y (~f = tt) = 
—i exp(— i^o~ y /2)\ 1=T: is applied which effectively changes 
the sign of the interaction in Eq. (fTJ). Hence the observ- 
able signal changes from s(t) in Eq. |4| to 

s n (t) = (UDl\n/2yrfo-£(t)RDl { (Tr/2)\V (7a) 
R = a c y f (S n t)a c y f (S n - 1 t)...o- c y i (S 2 t)<T c y f (5 1 t) . (7b) 

The evaluation of s n (t) is based on the same identities 
as the one of s(t) except that it is algebraically more 
involved. The result is cast in the form 

s n (t) = cos(2ip n (t)) exp(-2x,j(*)) with (8a) 

Vn(t) = J Q ^-X n (ut)dLU (8b) 

Xn(t) = J ^-coth{/3uo/2)\y n (Lut)\ 2 dLu , (8c) 
where the factor x n (z) in the integrand of the phase reads 

n 

x n (z) = (-1)" sin(z) + (-l) m+1 sin(z(5 m ) (9) 

m=l 

and the factor y n (z) in the integrand of Xnif) reads 

n 

y n (z) = 1 + (-l) n+ V z + 2 {-l) m e lzS ™ . (10) 

m=l 

The phase is less harmful since its influence on the 
signal is only quadratic in a. Hence we focus on Xn- 
The procedure discussed so far is based on n equidistant 
pulses @, 0, @ with S m = m/(n + 1) yielding 

\y n (z)\l q = 4tan 2 (z/(2n + 2))cos 2 (z/2) Vn even (11a) 
| y „(z)| 2 q = 4tan 2 (z/(2n + 2))sin 2 (z/2) Vn odd .(lib) 

For small values of z/(2n + 2) these functions rise 
quadratically like (l/2)z 2 /(n + l) 2 omitting rapid oscil- 
lations. Even without performing the integrations in (|5]l 
one can read off two features: (i) a large number n of 
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pulses is advantageous. The time scale tc is prolonged 
like tc — > (n+l)tc- (ii) no further suppression is achieved 
since the power law oc z 2 = (cot) 2 remains unchanged. 

Results of the equidistant 7r-pulse sequence are shown 
in Fig. [2] as dashed-dotted lines. Clearly, a shift to the 
right is discernible reflecting the growing factor n + 1. 
But no significant change of the slopes occurs. It is still 
important to have a weak coupling between qubit and 
bosonic bath to store the qubit for a significant time. 
For instance 100 pulses make it possible to store the qubit 
up to an error of 10~ 4 for » 5tc at a = 0.25 while for 
a = 0.001 it may be stored for w 60ic- 
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FIG. 2: Signal (|8]) vs. time for various numbers of pulses at 
T = 0. Solid lines for the optimized sequence, dashed-dotted 
lines for the equistant sequence (see main text). From top to 
bottom the curves refer to a = 0.25, 0.1, 0.01, 0.001. 

Now we pose the question whether the sequence of 
pulses can be optimized, for instance by exploiting the 
freedom of choosing the instants of the pulses. The po- 
tential of non-equidistant pulse sequences was recently 
demonstrated by concatenated pulse sequences [14J. In 
our work, we aim at finding the optimum pulse sequence 
with respect to canonical requirements. 

Inspecting (fTU|) one realizes that y n (0) = always and 
that there are n free parameters {5j}. So one may require 
that n additional conditions are fulfilled. We use this 
freedom to make the first n derivatives ?/n (z)| z =o with 
j G {1,2, ...n} vanish. Nicely, the resulting equations 
have a simple analytic solution for n 7r-pulscs 



Sj =sin 2 (7rj/(2n + 2)) 



(12) 



This is the main result of this Letter. The resulting factor 
in the integrand yields 



|ynO)lop = 



^2 (-l)- 7 'e (iz/2)cos(7rj/( ™ +1)) 



j——n—l 



16(n+l) 2 J 2 +1 (z/2) 



(13a) 
(13b) 



where the second line with the Bessel function J„+i rep- 
resents a very good approximation valid for z/(2n + 
2) < 1 up to exponential corrections. Note that [n + 
l)J 2 +1 (z/2) cx [z/(2n+ 2)] 2 "+ 2 manifesting the effect 
of the vanishing leading derivatives. From (fT3"|) follows 
that the integrand stays extremely small up to a certain 
value of z of the order of unity implying that decoher- 
ence hardly takes place up to a certain time t op given by 
t op ~ (n+l)ic- Beyond this time it sets in very abruptly. 

How does our finding compare to known results? For 
n = 2 we retrieve from Eq. |T2] <5i = 1/4 and 62 = 3/4. 
This means that our pulse sequence with the initial ir/2- 
pulse about and two 7r-pulses about a v reproduces the 
CPMG-cycle [6( which is widely considered for decoher- 
ence suppression 0, EH]. For all n > 2, Eq. [T2l predicts 
so far unexplored pulse sequences with a better potential 
for decoherence suppression. 

Fig. [2] depicts the features of the optimum sequence. 
Clearly, the lines are shifted to the right reflecting the 
factor (n + 1) in t op in parallel to the effect of equidis- 
tant pulses. In contrast to equidistant pulses the opti- 
mized pulses lead to steeper and steeper slopes on in- 
creasing n implying that the behavior for different cou- 
plings a becomes almost indistinguishable. This feature 
is extremely advantageous because it means practically 
that even a large coupling a does not harm a long stor- 
age time. For instance for 100 pulses the qubit may be 
stored for « 200ic independently of the value of a. The 
possible storage time is by a factor of 40 better than for 
the equidistant scheme for a = 0.25. For a = 0.001 the 
improvement is still about a factor of 4. Clearly, the 
improvement is most striking for large values of the cou- 
pling. Recurring to the estimate of tc = lps J3} we see 
that 100 pulses allow us to extend the storage time to 
about 200ps. 

Another way of looking at the optimized scheme Eq.[T2l 
is to ask how many pulses are needed to achieve a cer- 
tain storage time with an error below a certain threshold, 
say 10~ 4 . For a = 0.25, 5 or 6 pulses already imply a 
storage time of 5£d- For the same storage the equidis- 
tant scheme requires about 100 pulses. Keeping in mind 
that in practice each pulse will be imperfect it is certainly 
advantageous to work with a minimum number of pulses. 

Let us turn to temperature. In practice, no system will 
be at T = and in particular the favorable soft media 
will be operated at relative high T compared to the cutoff 
temperature Td = wd (setting fee = 1 = K). In our 
model, T = 1/(3 enters in the Eqs. (|6c|8c|) by the coth 
factor reflecting the thermal occupation of the bosonic 
modes. It deviates from its T = value of unity only for 
small frequencies. Small frequencies mean small values of 
z = uit so that the suppression of |2M(^)|op m this range, 
see Eq. (|13bj) . is particularly helpful. Finite temperature 
does not lead to any noticeable decoherence as long as the 
storage time is not too long, i.e. as long as it stays below 
t op . Indeed, curves at T ^ are indistinguishable from 
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the solid ones in Fig. [2] This holds already for a rather 
small number of pulses so that it is a highly relevant 
feature for experimental realizations. For the equidistant 
scheme the thermal effects are larger, in particular for 
large temperatures T Trj. The extreme insensitivity 
to thermal effects represents another essential advantage 
of the optimized scheme. 

The above derivations hold for arbitrary T, i.e. even 
for the classical limit T — > oo. Indeed, the same pulse 
sequences can be used for the classical Hamiltonian H = 
(f(t)/2)(j z where f(t) is controlled by Gaussian fluctu- 
ations determined by (/(£)) = and by (/(ii)/(*2)) — 
g(t\ — fa). The Fourier transform p(w) of g(t) is the 
power spectrum and p(u)/n replaces J(lu) coth(/3w/2) in 
Eqs. ()6c|8cj) while the phases <p and ip n are zero classi- 
cally. The other equations remain the same, in particular 
Eqs. (|12I13[) . This observation greatly increases the ap- 
plicability of our findings since many systems, not only 
bosonic baths, can be described for high temperatures by 
classical Gaussian fluctuations. 

The equidistant dynamic decoupling or iterated 
CPMG cycles have been realized experimentally, in par- 
ticular in NMR experiments. The detrimental influence 
of very slow nuclear spins on a solid-state qubit [l6| or 
on the electron spin in quantum dots [13] has been re- 
duced recently. A Rabi oscillation could be made vanish 
by realizing sequences of almost instantaneous 7r-pulses 
exploi ting the interplay between nuclear and electronic 



spins |18|. Krojanski and Suter demonstrated recently 



that even the decoherence of large quantum registers, 
realized by nuclear spins and their dipole-dipole interac- 
tion, can be significantly reduced by dynamic decoupling 
[l^ . To our knowledge, however, no optimized sequences 
obeying Eq. (|12p have been examined. 

An optimized sequence is by definition more efficient 
than a random one, cf. Ref. [2(|. But for large sym- 
metry groups it may be easier to use a random scheme 
than to optimize the pulse sequence. If more specific in- 
formation on the bath is available, cf. Ref. [15j, other, 
specifically adapted schemes might work more efficiently. 
Furthermore, we emphasize that hybrid techniques are 
attractive: Any other dynamic decoupling scheme may 
be improved by replacing the 7r-pulse or the CPMG cy- 
cle of two 7r-pulses by a suitable n > 2 sequence obeying 
Eq. HU The optimized design of real 7r-pulses of finite 
duration in the presence of bosonic baths, cf. Ref. [2l| 
for classical baths, is left for future research. 

In summary, we discussed strategies for suppressing 
the decoherence of physical quantum bits by dynamic 
decoupling. A promising way to optimize the sequence 
of 7r-pulses beyond the well-known CPMG sequence was 
analytically established. The comparison to equidistant 
pulse sequences revealed that the optimized scheme en- 
hances the possible storage time by up to almost two 
orders of magnitude. Alternatively, the number of pulses 



required to achieve a certain prolongation of the stor- 
age time can be much smaller (by a factor of 20 for 
strong coupling to the bosonic bath) than for the stan- 
dard equidistant scheme. Additionally, the optimized 
scheme is extremely insensitive to detrimental thermal 
fluctuations. So experimental investigations of the opti- 
mized scheme are called for. 
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